Fermi National Accelerator Laboratory 


FERMILAB-Pub-91/340-A 
November 1991 


Axino- induced baryogenesis 

//t /-72,-0/C 

Silvia Mollerach and Esteban Roulet a a 

Oj I 4- 


NASA/Fermilab Astrophysics Center 
Fermi National Accelerator Laboratory 
P.O. Box 500, Batavia, IL 60510-0500 


We consider the possibility that the baryon asymmetry is generated at low 
energies as a consequence of the axino decay. We analyse models in which the 
axino, the superpartner of the axion, is heavy and decays into gluinos at tem- 
peratures T ~ 1 GeV. If CP and B violating couplings for the quark superfields 
are included, the. subsequent decay of these out of equilibrium gluinos can act as 
seeds for baryogenesis. The required amount oi CP violation is well consistent 
with the bounds on the electric dipole moment of the neutron and the mechanism 
works even for low reheating temperatures after inflation (Trr ^ 10 4 GeV). 
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I. INTRODUCTION 


The standard explanation for the origin of the observed baryonic asymmetry 
77 = ns/n, ~ IO -10 requires B and CP violating processes occurring out of equi- 
librium during the early stages of the evolution of the universe. 1 These so-called 
Sakharov conditions are naturally met in some grand unified theories, which contain 
heavy Higgs (and vector) bosons with the desired couplings that may be left out of 
equilibrium by the fast expansion of the universe, producing the baryon asymmetry 
when they decay. 

The required conditions are in principle also attainable in the standard model, 
that allows for anomalous B violation, although in this case the out of equilibrium 
situation is difficult to realize at the low temperatures involved. Since the anomalous 
electroweak B violating (and B — L conserving) processes are in equilibrium up to 
temperatures T ~ even if they happen not to be relevant for the generation 

of 77, they are important because they may erase a baryon asymmetry generated 
previously (e.g. in a GUT scenario), unless B — L ^ 0 initially. 2 Furthermore, if the 
anomalous processes are combined with other sources of low energy L or B violation 
that are present in many extensions of the standard model (Majorana masses, R 
parity violating couplings, etc), they can eliminate any preexisting combination of B 
and L? 

Another possible problem for the GUT baryogenesis is that if the universe under- 
goes a period of inflation, as is required to solve the problems of flatness, horizon, 
monopoles, etc., the reheating temperature after the end of inflation may not be high 
enough to allow the heavy bosons to be thermally produced again (Trr ^ 10 12 GeV), 
in which case one is forced to devise an alternative mechanism leading to baryogenesis 
at low energies. In addition to the electroweak baryogenesis, 2,4 some other scenarios 
for this have been proposed. In particular, the supersymmetric extension of the stan- 
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dard model has been a playing ground for these attempts. Baryon violation can be 
obtained in these models by including the superpotential couplings gijk^idjdl, that are 
consistent with gauge and supersymmetry invariance and, if not combined with other 
L violating Yukawa couplings, are compatible with proton stability. Many new sources 
of CP violation are also usually present in supersymmetric models, 5 and the source of 
non-equilibrium has been searched in late-decaying particles. For example, in ref. [6] 
out of equilibrium squarks are produced by the decay of the inflaton itself and ref. [7] 
exploits the gauginos that result from the late gravitino decays. These out of equilib- 
rium particles can then generate the baryon asymmetry in their subsequent decay. In 
particular, the gravitino induced baryogenesis takes advantage of the possible over- 
abundance of gravitinos, that results from their very early decoupling at the Planck 
epoch and their survival until the decay time td ~ 8irMp/ml/ 2 ~ (10 TeV/m 3 / 2 ) 3 sec. 
To avoid disrupting the nucleosynthesis predictions of the light elements abundance, 
this scheme requires gravitino masses m 3 / 2 in excess of ~ 50 TeV. One possible in- 
convenience of this model is that if inflation dilutes the primordial gravitinos, their 
regeneration up to a sufficient amount, that is due to processes with gravitational 
strength, requires extremely large reheating temperatures, greater than ~ 10 15 GeV. 

In this paper we want to consider a new possibility that appears when one com- 
bines supersymmetry with the Peccei Quinn solution to the strong CP problem. 
In this frame, besides the axion the theory contains its superpartner, the axino a, 
which can have very important cosmological implications. 8 Since the axino couplings 
are suppressed by the PQ symmetry breaking scale F, that is constrained to be 
10 10 -j- 10 13 GeV, 9 axinos decouple at very early times surviving in large amounts. If 
axinos are stable, they would overdose the universe unless m 5 & 2 keV 8 and hence 
axinos heavier than this bound should be unstable. We will consider henceforth the 
supersymmetrized version of the invisible axion model of Kim 10 , in which axinos 
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are naturally heavy, with a mass that can be comparable to the gravitino mass. If 
m a ~ 1 TeV axinos will mainly decay into a gluon gluino pair with a lifetime 
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(This is similar to the 7 — ► <17 lifetime computed in ref. [12] but including the color 
factor.) For the allowed values of F these decays correspond to temperatures around 
1 GeV and thus do not interfere with nucleosynthesis at all, and one is tempted to 
use them as the source for the out of equilibrium gluinos that can act as seeds for 
baryogenesis. Furthermore, since axinos are regenerated much more efficiently than 
gravitinos, this mechanism works also for much smaller reheating temperatures than 
those required in the gravitino induced baryogenesis of ref. [10]. As it will also turn 
out, the CP violation necessary to account for sufficient baryon asymmetry generation 
can be obtained for a wide range of the parameter space in a way compatible with the 
experimental constraints such as the induced electric dipole moment of the neutron. 
One should also note that in the scenario to be discussed here the natural dark 
matter constituent is the axion itself, which is one of the preferred Cold Dark Matter 
candidates. 


II. AXINO COSMOLOGY 

We will consider the supersymmetric version of the heavy quark invisible axion 
model proposed by Kim. 10 In this model, a heavy color triplet Q is introduced which 
can be taken to be an 517(2) x Z7( 1) singlet. Its right-handed component Qr carries 
PQ charge +1, and the left-handed component Ql carries —1. The relevant Yukawa 
coupling is given by 

C Y = -/Ql^Qr ~ h.c. , (2) 

where a is the complex scalar field singlet of 517(3 ) x 517(2) x 17(1) in which phase 
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resides the axion field. The PQ symmetry is realized for a PQ charge —2 for <r. 

When the PQ symmetry is spontaneously broken at a scale F = (a), the quarks 
Q attain a mass /(<r). The main contribution to the axino mass comes from the one 
loop diagram represented in figure 1, that leads to 13 

ma = = 1 TeV (3) 

where A is a number of order one. For / ~ 0(1), as we will assume, the axino is 
naturally heavy and to avoid overdosing the universe it should be unstable. 

In this model, the axion mainly interacts with ordinary matter through a loop of 
heavy quarks Q coupling to two gluons. This vertex can be written as 11 

C.GG = « W (4) 

The corresponding supersymmetric vertex of the axino-gluon-gluino interaction is 
given by 
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If gluinos are lighter than the axino, this coupling allows for the decay a — ► gg with 
the lifetime given in eq. (l). 

In order to study the cosmological evolution of axinos, one should follow their 
thermal history. At high temperatures the axinos are kept in thermal equilibrium 
through the gluon mediated reactions qq «-► ag and qg «-► qa , which proceed at a rate 

*3 
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where n[T ) = 2£(3)T 3 /7r 2 is the number density of a relativistic specie at temper- 
ature T. When T becomes smaller than the expansion rate of the universe H, the 
interactions are no more effective to maintain axinos in chemical equilibrium and they 
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decouple. This happens at a temperature 


T d ~ 16tt 2 ~ 10 n Gevf — ) , 

15 V 90 7 M P a f \10 12 GeV/ ’ 
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where g, is the effective number of degrees of freedom, that for temperatures higher 
than the superpartner masses is g , = 930/4. 

At temperatures larger than Tp, the axino number density equals the photon one. 
After the decoupling both are equally diluted by the universe expansion; however 
successive ann i hilation of species increase the number density of photons leaving un- 
changed that of axinos. Thus, at a given temperature T < To, and before the axinos 
decay, the ratio of the number density of axinos to photons is given by 

9 .{T) 
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The picture may be changed if the universe underwent a period of inflation. The 
exponential expansion dilutes the initial number density of particles to negligible 
amounts. When the inflaton decays, reheating the universe, a new thermal bath of 
particles is created. If the reheating temperature Trh is larger than Tp>, the thermal 
bath will contain also axinos and the previous scheme holds. However, if the universe, 
reheats to. a temperature smaller than Td, the density of axinos will never reach the 
thermal equilibrium- one. In this case, we can compute the amount of regenerated 
axinos integrating the equation 


£ + 4*-r.(n 
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where R denotes the scale factor of the universe. The resulting ratio of the number 
density of axinos to photons for T < Trh is 

na = g.m ( g-(Tp) Trh 
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where we have introduced Tp> from eq. (7). 
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When axinos become non-relativistic (T < ma), the energy density p& = m a n & 
associated to the axinos decreases with temperature as T 3 , while that of relativistic 
particles as T*. Thus, the axinos contribution to the total energy increases as the 
universe expands and they may become dominant before they decay. 

The temperature of the universe just after the decay process occurs, T r , can be 
computed equating r to H~ l . As the axinos decay into relativistic particles, assuming 


that the reheating is instantaneous, we obtain using eq. (1) 

X /4 O. 3 /2 AJ 1 /* 
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(actually the universe does not reheat but just cools more slowly during the axino 
decay, but this treatment is satisfactory for our purposes.) 

If the baryogenesis results from the decay of the axinos, with each axino decay 
giving rise to e baryons, then the ratio of baryons to photons at present (To ~ 2.7° K ) 
is given by 

( i2) 

<7«(^r) n~( (T-) 

where g*(T 0 ) / g,(T r ) takes into account the increment of the number of photons with 
respect to that of baryons due to the annihilation of species after the axino decay. 
The baryon asymmetry can be written using eqs. (12), (l) and (11) as 

g,{T 0 ) y/2-rra, ( 87ry.(T r )-Y! /4 y/M P m a p a (r ) 


tj ~ e 


J.(T.) 8((3) 


90 


F Ptot{T r ) 
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The last factor measures the amount that the axinos contribute to the total energy 
when they decay. If they are dominant, it is of order one. This turns out to be 
the case when there is no inflation (or if the reheating temperature is larger than the 
decoupling one) and for values of the PQ scale F larger than ~ 10 11 GeV. For smaller 
values, the axino contributes 
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In the case that the reheating temperature after inflation is smaller than Td , the right 
hand side of eq. (14) should be multiplied by [g*{To)l <7 »(Tr #)] 3 / 2 Trh/Td. 


III. BARYON ASYMMETRY 

In order to generate a non vanishing e, we will consider the same mechanism that 
was proposed in ref. [7] to produce the baryon asymmetry out of gravitino decays. 
For this we start from a general softly broken supergravity lagrangian requiring lepton 
number conservation. This allows for the B violating superpotential couplings of the 
517(2) quark singlet superfields 


9ijk u id C jd c k , (15) 

with i,j, k being generation indices and for definiteness we will take only 0332 = 
<7332 /47T to be non vanishing ( stringent constraints apply to the couplings involving the 
first generation because of their effects on heavy nuclei decay and n n oscillations 6 ). 
The relevance of this term is twofold. First, it allows for CP violation in the gluino 
decays, because the interference between the second diagram of fig. 2, that involves 
the trilinear scalar interaction <7332 A77i 3 / 2 £ c 6 c .5 c , and the tree level term is a source of 
CP violation if lm(Am 3 / 2 m§) ^ 0. The two diagrams of fig. 2 will then generate 
more t t than t t pairs (or viceversa, depending on the CP violating phase) in the 
gluino decays. As was estimated in ref. [7], the resulting asymmetry AFj = r(<7 — ► 
tt) — T(g —> tt) is 


AT* 

r* 


**332 
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(16) 


with 

G(x) = «(l-x)(l- fT -^ J bi(x + i-l)). (17) 

The other important aspect of the couplings in eq. (15) is that they allow for the 
A B = ±1 decay of the (anti)squarks produced, in the gluino decay. Hence, each 
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gluino will give rise either to 3 quarks or 3 antiquarks, and the imbalance between 
these two channels, resulting from the CP violation, leads to a baryon asymmetry 
per axino decay e ~ AT§/Tg. 

In order that the generated baryon asymmetry is not erased by squark mediated 
processes such as qq — ► we will take the neutralinos x to be heavy (m x > 100 
GeV), so that these B violating processes are very suppressed after the axino decay. 
Note that in this model the neutralinos decay through x or x ~ ¥ 3q, and the 

lightest supersymmetric particle is not stable. 

Looking back to eq. (13), we see that if the axinos contribute significantly to the 
energy density at the time they decay (this holds for Trh Td), the value of e that 
is required to produce q ~ 10~ 10 is only e ~ 10“ 8 — 10 -7 . One can check that the 
required amount of CP violation does not induce at the same time an electric dipole 
moment for the neutron d n in conflict with observations. In fact, taking for instance 
<*332 = 0.01, mg — 300 GeV and m q - = 100 GeV, the electric dipole moment generated 
by the phase in Amg turns out to be, using the results of ref. [14], four orders of 
magnitude below the experimental upper bound |d n | £ 10~ 25 e cm. 15 This also allows 
to achieve the required values of q in scenarios with lower reheating temperatures, 
just by invoking larger values of e that can be easily obtained. For example, in the 
'extreme’ case a 332 — 0.1 and with and induced d n ~ 10 -25 e cm, Trh may be as low 
as 10 4 GeV for F ~ 10 n GeV. 

As a summary, we have considered supersymmetric invisible axion models in which 
axinos are heavy, as can happen in the Kim model. Assuming that the supersym- 
metric masses satisfy < mg < ms ~ 1 TeV we have shown that axinos decay at 
a temperature ~ 1 GeV, producing out of equilibrium gluinos. Introducing CP and 
B violating couplings these decays can act as seeds for baryogenesis. The observed 
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baryon asymmetry can be generated even with small amounts of CP violation and 
also for low reheating temperatures after inflation. 

It is a pleasure to acknowledge M. Turner for an enlightening discussion. This work 
was supported in part by the DOE and by NASA grant NAGW-2381 at Fermilab. 
The work of S. M. was supported by “Fondazione A. della Riccia”. 

Figure captions: 

Fig. 1: One-loop contribution to the axino mass. 

Fig. 2: Tree and one-loop diagrams contributing to the gluino decay. Their inter- 
ference is the source of CP violation leading to T(p — * tt) ^ T(g — *• tt). 
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